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Abstract 

K»" ] We have performed three dimensional nonlinear fluid simulations of electron fluid turbulence at 

C^^ ■ nanoscales in an unmagnetized warm dense plasma in which mode coupling between wave function 

and electrostatic potential associated with underlying electron plasma oscillations (EPOs) lead 
^D ' to nonlinear cascades in inertial range. While the wave function cascades towards smaller length 

oo ; 

^^ ■ scales, electrostatic potential follows an inverse cascade. We find from our simulations that quantum 



diffraction effect associated with a Bohm potential plays a critical role in determining the inertial 



H I range turbulent spectrum and the subsequent transport level exhibited by the 3D EPOs. 
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I. INTRODUCTION 



Studies of collective phenomena at nanoscales in dense matters are of great importance 

...,„„.,„,,*„.,.... ,,,»..,.,_„aaaaa^--. 

;ors [6j, nano-electromechanical systems [7[, quantum- diodes [8|, nanotubes and nanowires 
9|], quantum free electron lasers [lO|, as well as astrophysical bodies [ll|, [l2|] and intense 
laser-solid density plasma interaction experiments [13] for x-ray and 7-ray sources. In dense 
plasmas, the electrons are highly degenerate and quantum mechanical effects (e.g. electron 
tunneling arising from the finite width of the electron wave function) play an essential role 
at nanoscales. Since degenerate electrons follow the Fermi-Dirac statistics, there appear 
new electron equation of state and new forces involving the quantum Bohm potential [ij] 



and electron-1/2 spin effects 15j in dense quantum plasmas. It then turns out that due 



to the intrinsic nonlinearities associated with the Fermi pressure law and quantum forces, 

nn n 

there exists possibility of localizing electrostatic |16l . llTj and electromagnetic [18[ wave en- 



ergies at nanoscales in dense quantum plasmas. Here we report simulation studies of three- 
dimensional (3D) electron fluid turbulence at nanoscales in an unmagnetized warm dense 
)lasma. It is found that 3D nonlinearly interacting electron plasma oscillations (EPOs) 



191] in a dense quantum plasma exhibit nano-structures and associated energy spectra that 



are markedly different from those reported earlier for the 2D case [17|. Furthermore, we 
stress that the present 3D turbulence properties of electron plasma oscillations in our dense 
quantum plasma are significantly different from those in a classical plasma 20|, l2l|, |22| . In 



the latter, strong electron plasma wave turbulence has been studied by invoking parametric 
interactions [23| and by using either multi-dimensional cubic nonlinear Schrodinger equa- 



tion 



21 



22| or Zakharov equations 20], which are different from our 3D coupled nonlinear 



Schrodinger-Poisson equations in dense quantum plasmas. 

The surge for studying numerous nonlinear processes in dense quantum plasmas lies in 
a hope to transfer information through localized nano-structures one is able to create and 
sustain in plasmas. The present work dealing with 3D electron fluid turbulence shares a great 
deal of knowledge with classical fluid turbulence 24, 
superfiuid turbulence involving the Bose-Einstein condensates (BECs) 
gases. Both in fluids and plasmas as well as in BECs, one encounters the phenomena 
of inverse energy cascades in which energy transfer from small scales sustains large scale 



251 ]. plasma turbulence [26|, |27|, and 
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291 ] in ultracold 



circulations/structures in the flow, and results in a steady-state inertial range with power 



law scaling, as was originally predicted by Kolmogorov, Kraichnan and Iroshnik 2J, |25| 



while the dynamical equations depicting inverse cascades in fluids and plasmas are the 



Navier-Stokes and Charney-Hasegawa-Mima equations |26l.l27l| , the energy cascade scenario 



in BECS is described by the Gross-Pitaevskii equation |28l. |29| . In Sec II, we describe model 
equation governing the dynamical evolution of 3D dense fermi quantum plasma. We also 
present conservation laws admitted by the set of 3D equations. In Sec III, nonlinear 3D 
simulation results describing turbulence in such system are described. Mode structures and 
corresponding Kolmogorov-like spectra are also discussed. Turbulence transport is described 
in Sec IV, and flnally the conclusions are contained in Sec V. 



II. MODEL EQUATIONS 

In dense quantum plasmas, the Wigner-Poisson (WP) model has been used to derive a set 



of quantum hydrodynamic (QHD) equations 30| in the mean fleld approximation. The QHD 



equations include the continuity, momentum and Poisson equations. The quantum nature 
appears in the electron momentum equation through the pressure term, which requires 
knowledge of the Wigner distribution for a quantum mixture of electron wave functions, 
each characterized by an occupation probability satisfying the Pauli exclusion principle. 
The quantum part of the electron pressure is represented as a quantum force — V0b, where 
the Bohm potential is [l^ 0^ = — (^^/2?7iey^)V^y^. Here h is the Planck constant 
divided by 27t, rUe is the electron mass, and Ue is the electron number density. Deflning 



the effective wave function ijj = Jne{r,t) exp[iS{r,t)/h], where V5'(r, t) = meUe(r, t) and 
Ue(r, t) is the electron velocity, the electron momentum equation can be represented as 



an effective nonlinear Schrodinger (NLS) equation 3^, in which there appears a coupling 
between the wave function and the electrostatic potential associated with the EPOs. The 
electrostatic potential is determined from the Poisson equation. We thus have the coupled 
NLS and Poisson equations, which govern the dynamics of nonlinearly interacting EPOs is 
a warm dense quantum plasma. 

In this paper, we carry out simulations of 3D NLS and Poisson equations in order to 
understand the properties of 3D electron fluid turbulence (involving nano-structures and 
associated electron transport) in a warm dense plasma. We flnd that nonlinear couplings be- 



tween different scales EPOs are responsible for creating small-scale electron density clumps, 
while the electrostatic potential assumes large-scale structures. The total energy associ- 
ated with our 3D electron fluid turbulence at nanoscales processes a characteristic spectrum 
which is a non- Kolmogorov-like. 

For our 3D electron fluid turbulence studies, we use the NLS-Poisson equations [16|, 13 



iVz^— + ifV^^ + (/?^- |^|^/^^ = 0, (1) 

L/t' 

and 

VV = 1*1^-1, (2) 

which are valid at zero electron temperature for the Fermi-Dirac equilibrium distribution. 
In Eqs. (1) and (2) the wave function \1/ is normalized by ^/n^, the electrostatic potential 
if by ksTp/e, the time t by the electron plasma period ^p}, and the space r by the Fermi 
Debye radius Xd- We have introduced the notations Xd = {eokBTF/nQe^Y^"^ = Vp/^pe and 
yH = fiujpe/ \/2kBTp^ where eo is the electric permittivity, ks is the Boltzmann constant and 
the Fermi electron temperature ksTp = {Ti /2r7ie)(37r^)^''^nQ , e is magnitude of the electron 
charge, and Vtpe = {riQe^ /eorrieY/'^ is the unperturbed electron plasma frequency. The origin 
of the various terms in Eq.(l) is obvious. The first term is due to the electron inertia, the 
if -term is associated from the quantum diffraction effect involving the Bohm potential, (y9\l/ 
comes from the nonlinear coupling between the scalar potential (associated with the space 
charge electric field resulting form oscillations of the electrons around immobile ions) and the 
electron wave function, and the fourth term in the left-hand side of (1) is the contribution 
of the 3D electron pressure {pe = meVprv'J^ /hriQ ) for the Fermi plasma with a quantum 
statistical equation of state. 

Equations (1) and (2) admit a set of conserved quantities [l3|: the number of electrons 
N = J "^^dV, the electron momentum P = —i J \1/* V\l/(iy, the electron angular momentum 
L = -iJ^*rx V'^dV, and the total energy £ = J[-^*HV^^ + | V^| 72 + {3/5)\^\^^/'^]dV, 
where dV = dxdydz. In obtaining the total energy £, we used the relation dE/dt = 
iiJ(\l/V\I'* — \E'*V\1/), where the electric field E = —Vf. The conserved quantities are used 
to maintain the accuracy of the numerical integration of Eqs. (1) and (2). We note that lin- 
earizing the latter one obtains the EPO frequency cu = (fipg -|- k'^Vp + h^k^ /Ami] i "which 
exhibits the dispersive behavior of the EPOs. In the short wavelength regime characterized 



by k'^ ^ AmlVp/h'^, one notices that the dispersion associated with electron tunnehng effect 
dominates over that involving the quantum statistical electron pressure. 

III. NONLINEAR 3D SIMULATIONS OF QUANTUM PLASMAS 

We have developed 3D ffuid code to investigate nonlinear interactions between multi- 
scales EPOs described by (1) and (2). Our 3D ffuid code is based on Fourier expansion 



of the bases using a fully de-aliased pseudospectral numerical scheme 311]. The nonlinear 
de-convolution of Fourier modes is performed by computing the nonlinear triad interactions 
/(x, t)^(x, t) = I]k=k'-k" /(k', t)g(k", t)6(k' — k") at each time, which survive for only those 
coupled modes which satisfy the Fourier triad constraint k = k' — k". These nonlinear 
interactions in Fourier space follow from a Kolmogorov phenomenology of spectral energy 
transfer, and are mediated predominantly by the neighboring modes. Furthermore, such 
nonlinear interactions in the local spectral space conserve constants of motion of Eqs. (1) 
and (2), as presented above. The temporal integration is performed by 4th order Runge 
Kutta method. The spectral distribution for turbulent fluctuations is initialized isotropi- 
cally (no mean fields are assumed) with random phases and amplitudes in Fourier space. 
The evolution variables use periodic boundary conditions. The initial isotropic turbulent 
spectrum was chosen close to fc~^, with random phases in all three directions. The choice 
of such (or even a flatter than —2) spectrum treats the turbulent fluctuations on an equal 
footing and avoids any influence on the dynamical evolution that may be due to the initial 
spectral non-symmetry. Note, however, that the local as well as global mean flows may 
subsequently be generated by self-consistently excited nonlinear instabilities. Finally, the 
algorithm employed in our 3D fluid code ensures conservation of total energy and mean 
fluid density per unit time in the absence of charge exchange and external random forcing, 
and it is massively parallelized using Message Passing Interface (MPI) libraries to facilitate 
higher resolution. The conserved quantities, described in section II, are used to monitor the 
numerical accuracy of our 3D quantum fluid plasma code. Time step during the numerical 
integration varies between 10~^ and 10~^. The code uses usual 2/3 de-aliasing of the Fourier 
modes, such that the largest k is determined typically by k^ax = (2/3) (1/2) A^, where A^ is 
the total number of modes in each direction. The time step used in our 3D simulations not 
only preserves the conserved quantities, but it also minimizes the aliasing errors. 



(a) l^r, H=0.4 



(b) <i>, H=0.4 




(c) 1^1^, H=0.01 



(d) O, H=0.01 




FIG. 1: Fluctuations in the electron number density resulting from a steady state turbulence sim- 
ulations of our 3D electron quantum plasma. Forward cascades are responsible for the generation 
of relatively small-scale fluctuations in a decaying 3D electron quantum plasma as shown in (a). 
Large scale electron flow structures are present in the electrostatic potential, essentially resulting 
from a merging of smaller scale fluctuations as shown in (b). The latter is known as the inverse 
cascade process. H = 0.4 is used in (a) and (b). Figures (c) and (d) depict the electron number 
density and the electrostatic potential for H = 0.01. 

The localized initial turbulent spectral distribution, concentrated at the lower wavenum- 
bers, evolves in time following 3D nonlinear electron plasma wave interactions. Since the 
initial energy is localized in the large scale fluctuations, the latter drive turbulent processes 
through migration of energy towards relatively small scales Consequently, larger eddies trans- 
fer their energy to smaller ones through a forward cascade. During the forward cascade 
process, each Fourier mode in the inertial range spectrum obeys the vector triad constraints 
231 imposed by the vector relation k + p = q. These nonlinear interactions involve the 
neighboring Fourier components (k, p, q) that are excited in the local inertial range turbu- 
lence. We have performed a number of simulations to verify the consistency of our results 
in a strong turbulence regime. In our 3D simulations, we have explored two dense plasma 



systems that are characterized by different physical parameters, viz. dense warm plasmas 



in the next generation laser-based plasma compression (LBPC) schemes 32|, and the su- 



perdense astrophysical bodies 33|](e.g. interior of white dwarf stars). It is expected that in 
LBPC schemes, the electron number density may reach 10^^ m~^ and beyond. Hence, we 
have Qpe = 1.76 x 10^^ s-\ ksTp = 1.7 x 10-^*^ J, hQpe = 1.7 x 10"^^ J, and H = 1. The 
Fermi Debye radius Xo = 0.1 A. On the other hand, in the interior of white dwarf stars, 
we typically have uq ~ 10'^^ m^^ (such values are also common in dense neutron stars and 
supernovae), yielding Qpe = 5.64 x 10^^ s-\ keTp = 1.7 x 10"^^ J, hQpe = 5.64 x 10^^^ J, 
H ^ 0.3, and A^ = 0.025 A. The numerical solutions of Eqs. (1) and (2) for H = 0.4 and 
H = 0.025 (corresponding to uq = 10^^ m~'^ and no = 10^^ m~^, respectively) are displayed 
in Fig. (1), which are the electron number density (the left figure) and electrostatic (ES) 
potential distributions (the right figure) in the {x, y, z)-cube. 

One of most notable features of our simulations encompassing 3D electron fluid turbulence 
is that it exhibits a dual cascade phenomenon, as shown in Fig. 1. The electron density 
distribution in Fig 1 shows a tendency to generate smaller length-scale structures, while 
the ES potential cascades towards larger scales. The co-existence of small and larger scale 
structures in turbulence is a ubiquitous feature of various 2D and 3D turbulence systems. For 
example, in 2D hydrodynamic turbulence, the incompressible fluid admits two invariants, 
namely the energy and the mean squared vorticity. The two invariants, under the action of 
an external forcing, cascade simultaneously in turbulence, thereby leading to a dual cascade 
phenomena. On the other hand, 3D MHD turbulence exhibits forward cascades of energy 
and an inverse cascade of magnetic helicity. In these processes, the energy cascades towards 
smaller length-scales, while the magnetic helicity in MHD transfers spectral power towards 
larger length-scales. By contrast, the fluid vorticity in 3D hydrodynamics is prohibited 
from an inverse cascade. The randomly excited 3D Fourier modes nonetheless transfer 
the spectral energy by conserving the constants of motion in /c-space. In freely decaying 
quantum electron fluid turbulence reported here, the energy contained in the large-scale 
eddies is transferred to the smaller scales, leading to a statistically stationary inertial regime 
associated with the forward cascades of one of the invariants. Decaying turbulence often 
leads to the formation of coherent structures as turbulence relaxes, thus making nonlinear 
interactions rather inefficient when they are saturated. It is to be noted further that the 
long scale flow generation in our 3D simulations is observed to be directly proportional to 
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FIG. 2: Energy spectrum of 3D EPOs in a forward cascade regime. A Kolmogorov-like k~^^'^ 
has been observed for H = 0.4. The spectral index changes with respect to the parameter H. 
Numerical resolution is 128^. 

the parameter H. Intermittent flows are thus generated for a small value of H, while strong 
and large scale flows in the ES potential are formed when the magnitude of H is large (see, 
e.g. Fig. ([1])). The physical basis of such observation can be elucidated from the following 
arguments. The parameter H, which is the ratio between the energy density of the EPOs 
and the electron kinetic energy density of a warm dense quantum plasma, is associated with 
a diffraction-like term in Eq. (1) i.e. i7V^\t'. In this term, the negative imaginary part of the 
complex evolutionary variable \Ef essentially determines the rate of dissipation corresponding 
to the smaller scales. The smaller is H, more the dissipation is concentrated at the smaller 
scales and vice versa. For a moderately higher magnitude of the H parameter, there exists 
a strong tendency in EPO's to dissipate the smaller and intermittent turbulent eddies. It 
is therefore this H parameter which essentially characterizes electron flows at nanoscales in 
our 3D simulations. 

While the power spectrum for nonlinear EPOs exhibits an interesting feature in our 3D 
simulations, its scaling is not universal and is determined critically by the parameter H. 



x10 



15 




FIG. 3: Time evolution of effective electron diffusion coefficients measured for different values of 
H. Interestingly, smaller values of H corresponds to an effective low diffusion coefficient. The 
latter characterizes the presence of small scale turbulent eddies which results in suppression of the 
electron transport. 

For instance, we find a 3D Kolmogorov-like power spectrum k~^^'^ in some range of H 
values as shown in Fig. ([2]). The corresponding omnidirectional spectrum thus exhibits a 
^-5/3 gQaling. Spectral index nevertheless changes with H, as noted also in the study of 
2D fluid turbulence 171 . However, the spectral slope in the latter was found to be close to 
the Iroshnikov-Kraichnan (IK) power law 3J, |35j fc~^/^, rather than the usual Kolomogrov 



scaling 



36l | k ^/^. The origin of the differences in the observed spectral indices resides with 



the nonlinear character of the underlying warm dense plasmas, as nonlinear interactions in 
2D and 3D systems are governed typically by different nonlinear forces. The latter modify 
the spectral evolution of turbulent cascades to a significant degree. 

IV. ELECTRON TRANSPORT CAUSED BY TURBULENT FIELDS 



We finally study the electron diffusion coefficient in the presence of small and large scale 
turbulent EPOs. The effective electron diffusion coefficient produced by the momentum 



transfer can be calculated from -De// = /o°°(P(''7 1) ■ P(r, t + t'))dt', where P is the electron 
momentum and the angular bracket denotes spatial averages and the ensemble averages 
are normalized to unit mass. The effective diffusion coefficient -De//, resulting from 3D 
structures, essentially relates the diffusion processes associated with random translational 
motions of electrons in nonlinear fields of localized EPOs. It is remarkable to note that the 
electron transport can be effectively suppressed when the magnitude of the parameter H 
is decreased. This is shown in Fig. ([3]). The outcome of Fig. ([3]) is contrary to our 2D 



results 171], where the effective electron diffusion was lower when the field perturbations 
were Gaussian and increases later rapidly with the eventual formation of longer length-scale 
structures. Unlike the 2D case, the electron diffusion in 3D case is suppressed eventually 
because of the presence of small scale turbulent fluctuations. However, the steady state dif- 
fusion coefficient increases with the increase of H. This can be understood as follows. Higher 
the H valuer, the stronger is the small scale damping. When the small scale fluctuations 
are smeared out, the dynamical evolution of 3D quantum electron plasma wave turbulence 
is predominantly governed by the large scale flows, which consequently lead to a higher level 
of transport. This finding is further consistent with the fact that the parameter H dictates 
the characteristic evolution properties of the 3D EPOs, as described above. 

V. SUMMARY 

In conclusion, we have presented results from 3D nonlinear fluid simulations of the elec- 
tron fluid turbulence at nanoscales in an unmagnetized warm dense plasma. The mode 
couplings between the electron wave function and the electrostatic potential associated with 
the underlying electron plasma oscillations (EPOs) lead to the onset of nonlinear interac- 
tions and subsequent cascades in the inertial range. We find from our 3D simulations that 
the dispersion effect associated with the quantum Bohm potential plays a critical role in 
determining the inertial range turbulent spectrum and the subsequent transport level ex- 
hibited by the 3D EPOs. For instance, a Kolmogorov-like k~^^^^ 3D spectrum is observed 
for H = 0.4, whereas the spectrum flattens out for a smaller value of H. Correspondingly, 
the electron transport is higher for the higher H values. Finally, the wave function cascades 
towards smaller length scales, while the electrostatic potential follows an inverse cascade. 
We reiterate that the present investigation of 3D EPO turbulence is a necessary prerequisite 



for understanding the complex dynamical phenomena occurring at nanoscales in warm dense 
plasmas, such as those in laboratory and astrophysical settings. 
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